Spin Relaxation Anisotropy in Two-Dimensional Semiconductor Systems 
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Spin relaxation is investigated theoretically in two-dimensional systems. Various semiconductor 
structures of both n- and p-types are studied in detail. The most important spin relaxation mecha- 
nisms are considered. The spin relaxation times are calculated taking into account the contributions 
to the spin-orbit interaction due to both the bulk inversion asymmetry and the structure inversion 
asymmetry. It is shown that in-plane anisotropy of electron spin relaxation appears in III-V asym- 
metrical heterostructures. This anisotropy may be controlled by external parameters, and the spin 
I/") ' relaxation times differ by several orders of magnitude. 

CM 

I. INTRODUCTION 

i-Sh . Spin dynamics in semiconductors has been investigated for more than four decades. During this period, effects 
leading to creation and disappearance of equilibrium and non-equilibrium microscopic magnetic moments have been 
discovered. The non-equilibrium spin is known to appear in III-V and II- VI compounds upon absorption of circularly 
polarized light. This way to create non-equilibrium spins is named the method of Optical orientation of electrons and 
nuclei, which proved to be the most effective in investigating the optical and kinetic properties of bulk semiconductor 
C$ ■ samples and heterostructures. The problem of the loss of the average microscopic spin is very important in analysis 
of experimental data and device applications. 

In bulk semiconductors, mechanisms of electron and hole spin relaxation have been studied in ample detail both 
theoretically and experimentally . It is known that the main mechanism for electrons in crystals lacking inversion 
symmetry is, in a wide range of temperatures, the kinetic mechanism proposed by D'yakonov and Perel'. Carriers 
lose their spin orientation owing to precession in the effective magnetic field caused by those terms which are cubic in 
1— 1 1 wavevector. For holes in the valence band, where the spin-orbit interaction affects the energy spectrum to a greater 
extent, the loss of spin is mainly associated with spin-flip in each scattering act (Elliot-Yafet mechanism). However, 
the corresponding spin relaxation rate is of the same order as the momentum relaxation rate, which makes studying 
this kind of spin relaxation complicated. 
£f*^ ' In two-dimensional (2D) semiconductor structures, the relative importance of the above mechanisms changes dra- 
, matically because of the appearance of new forms of energy spectrum and spin-orbit interaction. The size-quantization 
CM ■ effect suppresses the spin-orbit interaction for holes and increases its strength for electrons. Therefore, it becomes 
possible to study spin relaxation in p-type structures, and in n-type samples it is fundamentally different from that 
in the bulk case. 

The spin dynamics processes in 2D systems may exhibit natural, 'intrinsic' anisotropy owing to presence of the 
growth axis. Analysis shows that in the conventional heterostructures this anisotropy is not large, being only 50%. 
However, the authors of the present study have shown theoretically that the in-plane anisotropy of spin relaxation 
appears in III-V asymmetrical heterostructures [Q. This anisotropy may be controlled by external parameters, and the 
spin relaxation times differ by several orders of magnitude Q . This opens up broad possibilities for spin engineering. 

We consider two most important carrier spin relaxation scenarios: the D'yakonov-Perel' and Elliot-Yafet mech- 
anisms. The Bir-Aronov-Pikus mechanism, associated with exchange electron-hole interaction, is essential only in 
p-doped samples, which are beyond the scope of this review, together with the excitonic spin relaxation. 
In general, spin relaxation equation can be written in the form 
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where Si is the spin density, and (l/iy ) is the component of a second rank tensor. The concrete form of the tensor 
(1/Ty) depends on the symmetry of the system under study. Since in a 2D structure the special direction exists 
(the growth axis) , spin relaxation processes are anisotropic. Symmetry of a single heterointerface of a semiconductor 
with zinc-blende lattice is C^v This leads to three different spin relaxation times. However if two interfaces of a 
quantum well are equivalent, then the structure symmetry increases up to D2d and in-plane spin relaxes isotropically. 
In semiconductors, spin relaxation processes are caused firstly not by external magnetic field, but due to peculiarities 
of a band structure and scattering mechanisms. These are processes which are considered in the present Review. 
We write the electron or hole Hamiltonian in the form 
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H = H o + V , 



(1.2) 



where Hq is the Hamiltonian of a system without weak spin-orbit interaction and scattering, and these two are 
described by small perturbation V. Let \nk) and E nk be the eigenstates and eigenvalues of Hq: 

(nk\H \n'k') = E nk 5 nn , <W . (1.3) 

Here n is the index of a 2D electron or hole subband, and k is the wavevector characterizing the motion in the 
heterostructure plane. In what follows, we assume that only the first subband of size quantization is populated 
and the energy dispersion is isotropic. Therefore, we omit the subband index and denote the energies as E k ■ It is 
convenient to extract the diagonal in the k part of the spin-orbit interaction: 

H so {k) = Vk k (1.4) 

(small spin-independent corrections from Vkk are assumed to be included into E k ). 

The dynamics of the electron spin density and that of the hole angular momentum density are determined by the 
time evolution of the respective density matrices. The diagonal in the k part of the density matrix, pk, obeys the 
equation 

-qT = r[#so(«),Pk] ~Y l*i [ 2 Vk'kPk'Vk'kl 6{E k -E k ,). (1.5) 
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This equation is valid for any type of spin-orbit interaction and scattering. Below we analyse it in the cases of electrons 
and holes separately. 

The review is organized as follows. Section II is concerned with 2D electron systems. Section III is devoted to 
p-type quantum wells. Conclusion discusses the possible ways to measure spin relaxation rates and to observe the 
spin relaxation anisotropy. 



II. ELECTRON SPIN RELAXATION 



A specific feature of electron systems is that the spin-orbit interaction is weaker than the momentum scattering. 
Therefore, it is convenient to represent the scattering amplitude as 

Vkk' = Tuw + V^' , (2-1) 

where X is 2 x 2 unit matrix and Tr(Vj^ k /) = 0. Since the spin-orbit interaction is weak, u 3> V . 

In the following, we assume a central elastic scattering. In this case, Ukk' is real and dependent on the electron 
energy E k and the absolute value of the scattering angle 9 — ip^i — y>k, where y?k is the angular co-ordinate of k in 
the heterostructure plane. 

Since iJgo(fc) and V^ k / are small perturbations, the spin relaxation times are much longer than the times in which 
the electron momentum distribution becomes isotropic. For this reason, it is convenient to represent the density 
matrix as a sum [Q : 

p = p + p', 7 = 0. 

Here the bar means averaging over ipy., and, hence, p depends only on E k . The anisotropic part of the density matrix 
is due to -ffso and V and, therefore, p' is sm all as compared with p- 



To second order in Hqq and V', equation (1.5) has the form 
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where g is the 2D density of electron states. It is noteworthy that the first two integrals in ( |2.2| ) are of first order in 
spin-orbit interaction, and the two last, of second order. 
The electron spin density is given by 

oo 

Se = g J dE k Tr(p^ . (2.3) 
o 

Therefore, to analyze the spin dynamics, we have to find p. Since Hso{k) and Vj^ k , are anisotropic in fe-space, the 
equations for p and p' are coupled and can be solved only for concrete interactions Hgo and V . 

In semiconductor heterostructures, spin relaxation can be caused by both H^o and V. The first scenario is known as 
the D'yakonov-Perel' mechanism, and the second, as the Elliot-Yafet mechanism. Even in the same heterostructure, 
the relative contributions of the mechanism may vary with such parameters as temperature or electron concentration. 
Since any of the spin relaxation mechanisms may be dominant, we investigate them separately. 



A. D'yakonov-Perel' mechanism 

Let us consider a system with spin-orbit interaction described by the Hamiltonian Hso{k). It is equivalent to a 
Zeeman term with effective magnetic field dependent on k. In the presence of scattering, the wavevector changes and, 
hence, the effective magnetic field changes too. Therefore, in the case of frequent scattering, electrons are subjected 
to a chaotically changing magnetic field. The spin dynamics in such a system has diffusion character, which leads to 
loss of specific spin orientation. This is the D'yakonov-Perel' spin relaxation mechanism p|, and it is the main spin 
relaxation mechanism in many III-V bulk semiconductors and heterostructures. It takes place even in the case of 
spin-independent scattering, when 

V = 0. 

In this case, in a time of the order of the transport time Tt r , the density matrix becomes isotropic, but the spin 
relaxation processes do not onset. One can obtain a kinetic equation for the spin density at times longer than Tt r [0: 
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Here 
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(2.5) 



are the Fourier-harmonics of the spin-orbit Hamiltonian, 

1 27r / d9 . .,. n . . „. 

— = Tfl f TT «kk' 2 (1 " cosnO) (2.6) 

are the transport scattering rates, and F±(Ek) are the distribution functions of particles with spin projection equal 
to ±1/2. 



Equation (2.4) is valid for 2D electrons with any spin-orbit interaction Hso(k). Let us consider an asymmetrical 
zincblende heterostructure. There are two contributions to Hso(k)- The first, the so-called bulk inversion asymmetry 
(BIA) term, is due to the lack of inversion symmetry in the bulk of the material of which the heterostructure is made. 
To calculate this term, one has to average the corresponding bulk expression over the size-quantized motion [ 5j . We 
investigate the heterostructure with growth direction [001] coinciding with the z-axis and assume that only the first 
electron subband is populated. The BIA term has the form 

H B i A (k) = 7 [(k 2 z )(<7 y ky - a x k x ) + a x k x kl - a y k y kl)} , (2.7) 

where we choose the x- and y-directions aligned with the principal axes in the heterostructure plane. Here (kf) is 
squared operator (—id/dz) averaged over the ground state, and 7 is the bulk spin-orbit interaction constant. It is 
seen that -ffeiA contains terms that are both linear and cubic in k. 
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In asymmetrical heterostructures, there is an additional contribution to the spin-orbit Ifemiltonian, which is absent 
in the bulk. It is caused by structure inversion asymmetry (SIA) and can be written as [ 



H S i A (k) = a (a x k y - a y k x ) 
where a is proportional to the electric field, E, acting on an electron: 

a = a>oeE . 



(2.8) 



(2.9) 



Here e is the elementary charge and ao is a second spin-orbit constant determined by both bulk spin-orbit interaction 
parameters and properties of heterointerfaces. It should be stressed that, in asymmetrical heterostructures, E is 
caused mainly by the difference of the wave function and band parameters at the interfaces, rather than by average 
electric field ||. 

i?siA contains terms linear in k, as also does -ffeiA- From equation (2.4) follows that the harmonics with the same 
n are coupled in the spin dynamics equations. This leads to interference of BIA and SIA terms linear in wavevector 
in the spin relaxation equations H. 

For Hso = -Hbia + ffsiA; the system has C2„-symmetry. Therefore, equations (|2.4|) can be rewritten as follows: 
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The times r z , r + and r_ are the relaxation times of a spin parallel to the axes [001], [110] and [110], respectively. 
If both the spin subsystems attain equilibrium before the onset of spin relaxation, then 



F±(E k ) = F (fi ± - E k ) , 



(2.11) 



where Fo is the Fermi-Dirac distribution function and fi± are the chemical potentials of the electron spin subsystems. 
If the spin splitting is small, i. e. 

- M-l < l/u+l) I A* — 1 

then the expressions for the spin relaxation rates 1/tj (i — z, +, — ) have the form 



1 

Ti 



J dE k (dF /dE k ) r,(fc) 
o 

CO 

/ dE k (dF /dE k ) 
o 



(2.12) 



where 
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(2.13) 
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Equations (|2.12|) , (|2.13D are valid for any electron energy distribution. If the electron gas is degenerate, then the 
spin relaxation times are given by 



- = r i (fc F ), 

n 

where fcp is the Fermi wavevector determined by the total 2D electron concentration ./V: 



(2.14) 



k F = V2ttN . 



(2.15) 

In this case, the scattering time T\ in equations ( [2.13] ) coincides with the transport relaxation time r tr , which can be 
determined from the electron mobility. 
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For nondegenerate electrons, the spin relaxation times are determined, in particular, by the energy dependences of 
the scattering times t± and T3. If Ti, T3 oc then t^Jt\ = const and 
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Here T is the electron temperature and &b is the Boltzmann constant. In the particular case of short-range scattering, 
v = 0, and ti = T3 are equal to r tr , which is independent of temperature. 

The spin relaxation tim es ar e very s ensitive to the relationship between two spin-orbit interaction strengths, 7(fc 2 ) 
and a. From equations ( [2.14 ), (2.16) follows that at low concentration or temperature, l/r z , 1/t_ and l/r+ are 
determined by respectively, the sum of squared 7(& 2 ) and a, their squared sum, and their squared difference. This 
may lead to a considerable difference between the three rates, i. e. to a total spin relaxation anisotropy if 7 (A: 2 ) and 
a are close in magnitude. In real III-V systems, the relations between Hbia and i/siA may vary, with Hbia or 
iJsiA dominant (l^Jll[l or these two terms comparable [l^ ]. In the last case, a giant spin relaxation anisotropy is 
predicted ||. 

The value of (fc 2 ) depends on the heteropotential and can be calculated for concrete asymmetrical heterostruc- 
tures [|| . The constant 7 is known for GaAs from optical orientation experiments [jy . Correct theoretical expressions 
for 7 and ao have been derived in terms of the three-band kp model |l3| , |l2| . The heterointerfaces make a contribution 
to ao in addition to that from the bulk JliJ]. At large wavevectors, ao starts to depend on k |ll|[l(|. Here we assume 
that the concentrations and temperatures are sufficiently low, allowing us to ignore this effect. 

In Ref. Q, the spin relaxation times were calculated for a III-V heterojunction and a triangular quantum well 
(QW) . The observance of spin relaxation anisotropy in all the three directions is predicted in a wide range of structure 
parameters and temperatures. 

It has been shown |p7j , |T8| , ^9| that inclusion of both the BIA and SIA terms, ( |2.7|) and (2.8), into Hgo leads to a 
spin-splitting anisotropy of the conduction band in fe-space in III-V semiconductor heterojunctions. However, the 
spin relaxation analysis performed in jll] ignored this effect. The authors of Ref. ]2(J] demonstrated that the BIA and 
SIA terms interfere in weak localization but are additive in spin relaxation. It was shown in Ref. that those terms 
in Hqo which are linear in the wavevector cancel out in spin relaxation as well. 

In a recent study the spin relaxation anisotropy was observed for uncommon (110) GaAs QWs. In this exper- 
iment, the spin relaxation in the growth direction was suppressed because of the 'built-in' anisotropy of the sample, 
resulting from the presence of heterointerfaces. By contrast, suppression of spin relaxation in the heterostructure 
plane was predicted in Ref. Moreover, all the three spin relaxation times are different in the last case, and this 
effect takes place in ordinary (001) heterostructures. 

We demonstrate that the terms in the spin-orbit Hamiltonian, which are linear in the wavevector interfere, and 
this leads to a huge anisotropy of the spin relaxation times. At high concentration or temperature, this effect starts 
to disappear owing to the predominance of the cubic in k terms in ifsoj which are only present in Hbia- However, 
the higher-order terms in Hqia are not forbidden by symmetry, either. These terms can also interfere with those in 
Hbia, and cause additional non-monotonic features in the dependences of the spin relaxation times on the structure 
parameters. 



B. Elliot- Yafet mechanism 



Let consider a 2D electron system without spin splitting of the spectrum. In this case, there is no Hso{k) term in 
the Hamiltonian, and the spin-orbit interaction occurs due to scattering only. 

In the heterostructures based on III-V semiconductors, the spin-flip scattering can be obtained in the Kane model 
taking into account a mixture of conduction- and valence-band states. The electron wave function has the form p2[ 

* s k(r) = c fc exp (ik ■ p) [u s (z)S + v sk (z) • R] . (2.17) 

Here s = +, — enumerates two spin states at a given k; S and R = (X, Y, Z) are, respectively, s- and p-like Bloch 
functions, and Ck is the normalization factor. The functions u s and t) s k are eight envelopes corresponding to the 
conduction and valence band, respectively, related by 

Vsk(z) = i (AK - iBa x K)u s (z) . (2.18) 
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Here we introduce a 3D vector K — (fc, k z ) with k z = —id/dz, and two constants: 



A = ^ ( J- 
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(2.19) 



where £J g and A are the energy gaps between the bands Tq and Ts, and T$ and IV, respectively, and P is the Kane 
matrix element. The envelope function for the conduction band is 



u s (z) = ip(z) w s 



(2.20) 



where ip(z) is the wave function of size-quantization, and w + =f , w_ =| are the spin functions. The normalization 
factor is given to second order in AK and BK by 



c k = l-(A 2 +2B 2 )(k 2 + (k 2 z )), 



(2.21) 



where (k 2 z ) = J dz(dip/dz) 2 . 

The sc atter ing matrix element (sk\V\s'k r ) is a 2 x 2 matrix with respect to the indices s and s' . Therefore, in the 
notation (|2.l[) we obtain 



uo(k — fc')(fc 2 ) — I dzuo(k — k' , z)(k z (pY 



w k k' = u Q {k - k')[l + (A 2 + 2B 2 ){k ■ k' - fc 2 )] + {A 2 + 2B 2 
V£ k ' = (2^5 + B 2 ) |y dzw (fc - fc', z) <p ifc z (^ [<r x (fc + k% + iu (k - k') a ■ (k x fc') j 



(2.22) 



Here 



wo(<Z,z) = J dpV{p,z)exp (iq-p), 
with y(p, z) being the scattering potential, and 

"0(9) = dz uo(q, z) p 2 (z) . 



(2.23) 



(2.24) 



The difference between Ukk' and uo(fc — k') is due to the spin-orbit interaction, similarly to the 3D case |23| ]. 

The two terms in V' are different: the part containing <j x and a y is linear in fc, but the term with a z is ~ Ar. Hence, 
at fc 2 <C (fc 2 ), the first part is dominant, and V^ k , can be represented as 



where 



V4' = v{k - k') [ax{k + k% , 



v(q) — (2AB + B 2 ) J dzu a (q, z) ip ik z (p . 



(2.25) 



(2.26) 



It is noteworthy that a nonzero contribution at k — k' is the impurity-induced SIA term physically meaning the force 
exerted on electrons by scatterers. This term is to be added to -ffso, and it manifests itself in the D'yakonov-Perel' 
mechanism only. 

It is seen from equation ( 2.25 ) that 

^kk< = ^k<k 

for central scattering, whereas v(k — k') = v(Ek, — VkD- Therefore, averaging of Eq. fl2.2| ) over ip^ gives 
2tt 
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n ^kk' P ^kk' - u kk' KPk' ~ Pk) V kk> + ^kk' (Pk' - Pk)\ 



(2.27) 



It can be shown that the term in square brackets is proportional to the integral 

d0u{e)v{9) sin (9 
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which is zero for central scattering. 
Using Eq. (2.3), we obtain 
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It can be seen that the Elliot-Yafet mechanism results in a 50 % anisotropy of spin relaxation times. 
To second order in A/ E s , we have for a Boltzmann gas: 
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(2.30) 



Here the wavevector transfer in the case of elastic scattering is described by 



Q = 



8mE k 



sin — 
2 



and account was taken of the fact that, at a given accuracy. 

Ti 2 



p2 
E a 



The factor r is on the order of Ei/E g , where E\ is the energy of the first level of size-quantization. For an infinitely 
deep rectangular QW of width 100 A, E g = 1.5 eV, and for randomly distributed <5-scatterers, r — 1.2 x 10~ 2 . Hence, 
for the parameters of GaAs (A = 0.34 eV, m = 0.067mo) and T — 77 K, we have T tI /T zz ~ 10 -5 . 

In addition to the scattering processes leading to spin relaxation considered above, there is a short-range interaction 
making possible scattering of electrons from the T$ or Tq band to I\ band H . The same mechanism takes place in 
the 2D systems as well p3]. 



III. HOLE SPIN RELAXATION 



A specific feature of p-type QWs is a strong spin-orbit interaction: the states at the top of the valence size- 
quantized subbands have a certain projection of the angular momentum on the growth axis. Therefore, ordinary 
spin-independent scattering leads to hole angular momentum relaxation even in a symmetrical QW, when there are 
no additional spin-orbit terms: i?so = 0. In other words, the Elliot-Yafet mechanism is very important for holes. 

In a symmetrical p-type QW, each hole level is doubly degenerate. The wave functions have the form 



* s k = exp (ik ■ p) F s (k, z) , 



(3-1) 



where s — 1, 2 enumerates the degenerate states, and Oz is the growth direction. 

In a QW based on a III-V semiconductor, the hole wave function is a superposition of those four states of the 
valence band top, which correspond to the angular momentum projections on Oz equal to 3/2, 1/2, —1/2 and —3/2. 
In the basis of these states, F s (k, z) are four-component columns which can be represented as p5| 



-vo C(z) 
ivi S(z) e^ k 
-V 2 C{z) e 2lVk 
iv 3 S(z) e 3 ^" 



iv 3 S(z) e- 3 ^ 
v 2 C(z) e~ 2lVk 
iv x S(z) e-^" 
vo C(z) 



(3.2) 
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Here C(z) and S(z) are, respectively, symmetrical and asymmetrical functions of the coordinate z, and the real 
coefficients v m (m = ... 3) for the ground subband are given by 

v m ~ k m , k^O. (3.3) 

An analog of the electron spin density in experiments on circular polarization of the luminescence in low magnetic 
fields is the pseudospin of holes, Sh- Its value is proportional to the degree of circular polarization of radiation at 
relatively weak hole orientation. The hole pseudospin density is given by 

S h = Y / Tr[p u J(k)}, (3.4) 

k 

where J(k) is the pseudospin operator for the first hole level of size-quantization. In the basis (3.2), J l (fe) has a form 
of 2 x 2 matrix with the elements 

4,,(fe) = J dzF](k,z)J i F s ,(k,z), (3.5) 

where J % are 4x4 matrices in the basis of the four states of the valence band top |2y] . Calculation yields 

J z (E k )=a zJz , (3.6) 



where 



J x +ijy ( ae- 2 ^" 



j x - u y ( o & e - 4 ^<= 

2 ~ V a ^ iiPh 



3 n 1 rj \ f (j . . ( 1 n 3 



(3.7) 



(3.8) 



ME k ) = ^v 2 Q - -v 2 2 j J dzC 2 (z) + i-vf - -vlj J dzS 2 {z) , (3.9) 
i(E k ) = v\ J dzS 2 {z) - V3v a v 2 J dzC 2 (z) , (3.10) 



b(E k ) = V3v lV3 / dzS 2 (z) -v 2 dzC 2 {z) . (3.11) 



It can be seen that J z depends on the energy E k , whereas J x and J y depend on both E k and the direction of k. 
Therefore, we have to find not only the isotropic part of the density matrix, p, but also the 2nd and 4th Fourier- 
harmonics. 

J 1 have complicated energy dependences. For instance, j z (E k = 0) = 3/2 for the first heavy-hole subband of 
size-quantization. The dependence j z on k is given in p8| for different size-quantization subbands. It follows from 
equation ( [3.3| ) that a ~ k 2 and b ~ k 4 for small k, i.e. J x , J y are nonzero owing to nonparabolicity effects only. 

A. Elliot- Yafet mechanism 

Let us consider scattering from randomly distributed short-range impurities 

V(r) = J2Vo5(r-Ri). (3.12) 



Calculating the scattering amplitude V ss >(k,k r ) and using Eq. (1.5), we obtain the following equation for the hole 
density matrix 
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(3.13) 



Here, the angular brackets denote averaging over Ri, and account is taken of the fact that, for (3.2) and (3.12), 

(V u ■ V 12 ) = (V22 ■ V 12 ) = 

and I V ss ' I depends on 9 = ipu — ipk>. Note that the spin-orbit interaction in p-type QWs is strong and, therefore, the 
division (2.1) makes no sense, and we use the total scattering amplitude Vkk'- 
We expand the density matrix in a series 



71— — OG 



The spin-independent part obeys the equation 



dt 



2tt 



-g(E k 



d6 
2^ 



(|Vii| 2 + |^i 2 | 2 )(l-cosn0)/ n 



(3.14) 



(3.15) 



It is seen that the total number of particles with a given energy fo(E k ) is conserved, which is correct for the elastic 
scattering under study. The Fourier-harmonics /„ with n > 1 relax with the corresponding transport times. 
For the z-component of ae„ we have 



dt 



~9(E k )j ^ [(|y 11 | 2 )(l-cosn0) + (|l/ 12 | 2 )(l+cos^)] 



(3.16) 



It is important that for small k, 



$V U » 9fVii » \V 12 \ 



Therefore, the hrst term in square brackets is close to the inverse transport time, and the second is much smaller and 
describes spin relaxation. It is seen that the average pseudo-spin with energy E k , se§, relaxes with time much longer 
that r tr . However, higher Fourier-harmonics change much faster. 

The same is true for ee x ' v . This can be readily shown for small k, when we can neglect V12 compared with Vu. In 
this limit 



d{sf n ± jeep 
dt 



fan i i&n) ' 



(3.17) 



2tt 

T 



/d6 
^ [{\Vu I 2 ) (1- cos n9) + 2((%V 11 ) 2 ) cos ri^2(»%. %Vn) sin nB] . (3.18) 



This means that all harmonics relax with time ~ Tt r except that with n = 0. The zero harmonic determines the hole 
Cooperon, important in the weak localization effect. Its relaxation is caused by the small term | V12 1 p6| , p7[ : Ho wever, 
the in-plane pseudospin components J x , J y are determined by the non-zero harmonics ee%'% [see equations (3.7), (fDj})] 
and, hence, their relaxation is fast. Therefore, one can investigate only the dynamics of Sh, z in optical experiments. 
The corresponding equation is 

00 

JdE k g*(E k )[F + {E k )-F_{E k )]j z {E k )§ § (\V 12 \ 2 ) 

-qT--j- 55 Sh ^- ^- 19 > 

/ dE k g(E k ) [F+(E k ) - F-(E k )] j z {E k ) 


In Ref. ps}] , the value {2tt /h)g(E k ) § d6 /2-k{\Vi 2 \ 2 ) was calculated as a function of k for different scattering potentials. 
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B. D'yakonov-Perel' mechanism 



If we take into account the splitting between the pseudospin states in the valence subbands, which may be caused 
by both BIA and SIA (in asymmetrical heterostructures), then the operator Hso f° r boles. For the ground 
subband, in the basis of the two degenerate states, it is a 2 x 2 matrix. Therefore, the problem is equivalent to that for 
electro ns, a nd we can obtain the equation for the density matrix dynamics. Expanding th e den sity matrix in Fourier 



series ( 3.14 ), we obtain expressions for /„ and ae„. The former coincides with equation ( 3.15 ), and the latter is as 
follows 

= ~^9{E k ) j> — r„Tr([tf_ nj [H n , a^ai) < . (3.20) 



According to equations (3.6)-(3.£), in order to investigate the hole pseudospin dynamics, we have to average over E k 
the equations with n — 0, ±2 and ±4. It can be seen that, similarly to the Elliot-Yafet mechanism, relaxation of Sh,x 
and Sh, y is fast, taking a time ~ rt r . This is due to the anisotropic nature of the operators J x and J v . The only 
slowly relaxing component is Sh, z - The corresponding equation has the form 

oo 

„„ . oo / dE k g(E k ) [F+(E k ) - F_(E k )] Jz (E k ) Tn Tr{[H_ n) [H n , a z ]]a z ) 

lf=-iE 5 °° S h , z . (3.21) 

2h »=— JdE k g(E k )[F + (E k )-F_(E k )]j z (E k ) 
o 

In [29| ], correct equations were obtained for relaxation of the zero harmonics of the density matrix (for /o and 
ae ). Therefore, the slow time calculated there describes the relaxation of Sh, z only. The relaxation time of aeg' y was 
incorrectly associated with that of Sh,x and S^.y 

To calculate the hole pseudospin relaxation time, a knowledge is necessary of the concrete form of i?so for a 
structure studied. The BIA contributions can be obtained by averaging the corresponding bulk terms p2| . Moreover, 
the SIA contributions exis t in a symmetrical heterostructures |3G] . All of them lead to hole pseudospin relaxation in 



accordance with equation ( 3.21 ) 



IV. CONCLUSION 

The predicted spin relaxation anisotropy can be observed in time-resolved measurements similar to those in Ref. [ pH . 
In steady-state experiments, the in-plane electron spin relaxation anisotropy can be investigated by means of the 
Hanle effect. To obtain the spin relaxation times, account should be taken of the fact that the Lande g-factor has 
not only diagonal in-plane components (g X x), but also off-diagonal ones (g X y) [0- The degree of photoluminescence 
polarization in a magnetic field B _L z is described by the following expression 

P{B) = , (4.1) 

1 + [^B (gxx ± 9xy) B/H} 2 T z T T ' 

where the upper and lower signs correspond, respectively, to the experimental arrangements with £?|| [110] and B\ | [110] 
{lis is the Bohr magneton). By changing the structure asymmetry (e. g., by applying a gate voltage), one can modify 
the Hanle curves owing to the in-plane spin relaxation anisotropy. 
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